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Ov| . Abstract 

^^ . The final data of the Gallex and GNO radiochemical experiments have been searched 

QQ I for possible time modulations of the solar neutrino capture rate using the Lomb- 

^^ ■ Scargle periodogram and the maximum likelihood methods. Both analyses, per- 

formed using the raw data, do not support the presence of a time variability with 
f^ I characteristic periods resembling those of rotation of the solar magnetic fields, as 

^T ' previously suggested in the literature by some authors. In this context, the poten- 

tial sensitivity of the Lomb-Scargle method to radiochemical ^^Ga data has also 
been explored with simulated data. This allowed to set an exclusion plot in the 
frequency/amplitude plane from the Gallex-GNO dataset. 
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1 Introduction 



In recent years it has been suggested [1,2,3,4,5,6,7] that solar neutrino flux 
measured by CI and Ga radiochemical experiments and by Super-Kamiokande 
shows a time variability with characteristic periods resembling those of rota- 
tion of the solar magnetic flelds. Since this could be a very strong hint for 
sub-leading Resonant-Spin-FIavour-Precession (RSFP) of pp neutrinos occur- 
ring in the solar convection zone [8] and for a non-negligible neutrino transi- 
tion magnetic moment, it is of paramount importance to perform an extensive 
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search for possible time modulations in the '^^Ga data using different methods 
and exploiting the total information of the available solar runs. 
The Gallex [9] and GNO [10,11,12] radiochemical experiments at Laboratori 
Nazionali del Gran Sasso, Italy, monitored the low energy solar neutrino flux 
using a 30-ton gallium detector via the inverse f3 reaction ^^Ga(z/e,e~)^^Ge. 
GNO is the successor project of Gallex, which continuously took data be- 
tween 1992 and 1997, and has been in operation from April 1998 to April 
2003. The gallium detector at Gran Sasso consists of 30.3 tons of natural gal- 
lium in the form of 103 tons of gallium chloride (GaCls) acidic solution. The 
energy threshold of the ^^Ga(z/e)e^)^^Ge reaction, 233 keV, is sufficiently low 
to make gallium detectors sensitive to the fundamental pp neutrinos. Presently 
only radiochemical experiments based on ^^Ga are sensitive to these neutrinos. 
^^Ge nuclei produced by z/g interactions are unstable and transform back to 
^^Ga by electron capture, with a life time r of 16.49 days. The gallium target 
is exposed to solar neutrinos for a time texp = iEOE-'^soE (SOE = start of ex- 
posure, EOE = end of exposure) which is typically 3 or 4 weeks. After this, the 
accumulated ^^Ge (~8 atoms for a 4- week-exposure) is extracted from the tar- 
get solution and converted in germane gas GeH4. The germane is added with 
low-activity xenon and is used as counting gas for a miniaturized proportional 
counter, which is able to detect the electron capture decays of ^^Ge. The inter- 
action rate of solar neutrinos with the target is measured, in the hypothesis it 
is constant in time, by counting the decay signals of ^^Ge in the proportional 
counter. For details about the experimental procedure and data analysis see 
Ref. [10]. The final result from Gallex (65 solar runs, 1594 days of exposure 
time) is 77.5±6.2±4.5 SNU. ^ From April 1998 to April 2003, 58 solar runs, 
with exposure time of 3 or 4 weeks, were successfully performed within the 
GNO project, for a overall time of 1713 days [10,11,12]: the resulting solar neu- 
trino interaction rate is 62.9l5;|±2.5 SNU. The combined Gallex/GNO rate 
(123 solar runs, 3307 days of exposure time) is 69.3±4.1±3.6 SNU. Results 
from the single solar runs are displayed in Fig. 1. 

In this paper, possible significant time periodicities in the solar neutrino data 
are looked for using the Lomb-Scargle periodogram method (which has already 
been used in the literature for this kind of analysis); the potential sensitivity is 
also tested with simulated data sets. The search for time periodicities is then 
repeated by applying the maximum likelihood method to the candidate ^^Ge 
decays in the single solar runs; the constant solar neutrino rate hypothesis is 
tested and a likelihood spectrum is derived. 



1 SNU = 10 ^^ interactions per second and per target nucleus. 
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Fig. 1. Measured solar neutrino capture rate (SNU) in the 123 solar runs of 
Gallex/GNO [12]. 

2 Lomb-Scargle periodogram method 



2.1 Analysis of the solar runs 



The search for periodicities in the Gallex/GNO solar neutrino data can be 
performed using the Lomb-Scargle periodogram method (see Ref. [2,6,13]). 
Given a set of A^ data points hi = h{ti) at the times {ti}, i = 1 . . . N, with 
mean and variance h and a^ respectively, the Lomb-Scargle periodogram is 
defined as [14,15] 
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where uo = 2ttv is the angular speed and v is the frequency; the offset 9 is 
defined by the relation 
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and makes the periodogram -2(z/), which is also called power, independent of 
shifting all the times {ti} by any constant [15]. 
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Fig. 2. Lomb-Scargle periodogram for the 84 solar runs of Gallex and GNO I [9,10] 
in the frequency range 10 — 20y~^. This plot corresponds to Fig. 2 of Ref. [2] 

Following the procedure used in Ref. [13], the periodogram for the Gallex/GNO 
data presented in Ref. [2] is first reproduced. Fig. 2 shows the Lomb-Scargle 
periodogram obtained from the analysis of the first 84 Gallex/GNO solar runs 
[9,10], that can be directly compared with Fig. 2 at the page 1369 of Ref. [2]. 
The most relevant features of the original periodogram are correctly repro- 
duced: the peaks (and in particular the main one, at frequency ~ 13.7 y~^) 
are approximately in the same position and the vertical scale of power, which 
fixes the statistical significance of the peaks themselves, is exactly equal. This 
level of agreement is completely satisfactory for the present testing purposes ^ . 
After this test of the Lomb-Scargle code, the same analysis has been extended 
to the whole set of Gallex and GNO data [9,10,11,12]. As in the preliminary 
test of Fig. 2, the results {hi} of the A^ = 123 solar runs are conventionally 
referred to the end-of-exposure {t^Q^;} times[2]. In order to test the region 
of major astrophysical interest, the frequency range taken into account spans 
from 
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where z/f. 
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5.19 y is the average Nyquist frequency. In the pe- 



2{tjv-ti) 

riodogram analysis all the data points {hj} are used with equal statistical 
weights and are corrected to take into account the expected l/(f' geometrical 
modulation of the neutrino flux, which is due to the variation of the Sun-Earth 
distance during the year. If Zmax is the highest peak found in the periodogram, 
one has to extablish if this is compatible with statistical fluctuations, i.e. one 



^ Minor differences in the positions or in the powers of the single peaks, that are 
probably related to different used approximations, are not significant in this context. 



600 t 
500 
400 F 
300 r 
200 E- 
100 




J I I L 








10 15 

Maximum power 



o 
o 
^ 0,8 h 



^ 0.6 

£ 0,4 

0.2 





_] I I L. 



_] I I L. 



-' ' I L. 



10 15 

Maximum power 



Fig. 3. Distribution of maximum power z^ax (left) and false-alarm probability 
(right) for the null hypothesis (i.e. constant signal with random noise), assuming 
N = 123 data with times {tj} as in the real solar runs and frequency in the range 

has to test the null hypothesis Hq of a constant neutrino signal with random 
noise. The probability P{> Zmax\Ho) that statistical fluctuations give a peak 
in the periodogram larger than z^ax is the false-alarm probability of the null 
hypothesis and a small value of P{> Zmax\Ho) indicates an highly significant 
periodic signal. Following the procedure of Ref. [16], the false-alarm proba- 
bility has been estimated by generating 100,000 Monte Carlo data samples 
{hi}, i = 1 . . . 123, in the assumption of the hypothesis Hq. As the probability 
depends also on the detailed time spacing of the data, the {hi} have been 
referred to the same times {tj} as the real runs. The Lomb-Scargle method is 
then applied to each data sample, to get the spectrum z^u) and the maximum 
power Zmax- The resulting maximum power distribution and the false-alarm 
probability are displayed in Fig. 3; the Zmax power values corresponding to 
50%, 5%, 1% and 0.1% false-alarm probabilities are respectively 6.69, 9.12, 
10.57 and 12.55.^ 

The Lomb-Scargle periodogram obtained from the analysis of the 123 solar 
runs of Gallex and GNO is shown in Fig. 4. The highest peak in the interval 
[^min,^max] IS fouud at Umax = 6.37 ± 0.04 y~^; the corresponding power is 
Zmax = 5.26 and the false-alarm probability for the null hypothesis is ~95%. 
The Lomb-Scargle periodogram for Gallex/GNO solar neutrino data is hence 
statistically consistent with the expectation of a constant interaction rate and 
it does not contain any hint of time modulation. If this same analysis is re- 
peated for the 58 solar runs of GNO only, the highest peak is found at the 
frequency Umax = 2.9 ±0.1 y^^, with false-alarm probability ~55%. 



^ The 50% false-alarm level, Zmaa;=6.69, corresponds to the average value of z„ 
expected from a null hypothesis signal. 
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Fig. 4. Lomb-Scargle periodogram for the 123 solar runs of Gallex/GNO in the 
frequency range [I'mimi^max]', dashed hues are the 50%, 5%, 1% and 0.1% false-alarm 
levels calculated from Fig. 3. 

2.2 Test of the sensitivity with simulated data 



In order to test the actual sensitivity of this kind of data to possible time mod- 
ulations, Monte Carlo repetitions {hi} of the Gallex/GNO experiments have 
been generated, in the hypothesis that the ^^Ge neutrino-induced production 
rate is modulated in time as 



n{t) = no[l + Asm{2nut + ip)], 



(4) 



with < A < 1. The mean production rate no is assumed to be Uq = 0.63 
atoms/day = 70 SNU, close to the experimental value. The statistical fluctu- 
ations on the single runs can be approximated with a Gaussian having rms 
a =~40 SNU. ^ The search for periodicities in the simulated data has been 
then performed using the Lomb-Scargle periodogram method in the frequency 
range [I'min, t^max] defined above. Being the Lomb-Scargle periodogram invari- 
ant for time translations [15], it can be assumed without loss of generality that 
ip = 0, since this is simply equivalent to shift the time origin. For each data 
set, the maximum power Zmax{A, v) was evaluated. Fig. 5 shows the Zmax{A, u) 
for frequencies ranging in 1-20 j~^ and for amplitudes A = 1, | and |. It can 
be seen from Fig. 5 that the method clearly indentifies the time modulation 
if the frequency is low and the amplitude is large enough. However, even for 
^4=100%, the sensitivity rapidly decreases for frequencies v higher than ~ 12 



^ The actual distribution of the results of single runs is not exactly Gaussian, 
because of the asymmetry introduced by Poissonian effects. However, it has been 
verified that the assumption of Gaussian fluctuations does not affect the final result 
in this particular application. 
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Fig. 5. Maximum power Zmaxi^,'^) vs. frequency obtained from the Lomb-Scargle 
analysis of simulated data, in the hypothesis time modulation with frequency u and 
amplitude A = 1, ^ and j. The statistical error on the single run is assumed to be 
40 SNU. Dashed lines are the 50%, 5%, 1% and 0.1% false-alarm levels. 

y^^: in this case, the modulation signal cannot be statistically distinguished 
from a constant signal with Gaussian fluctuation, i.e. it is impossible to reject 
the null hypothesis using this kind of data and analysis. 



2.3 Exclusion plot 



Since no positive modulation signal has been found with the periodogram 
analysis of the Gallex/GNO dataset in the region of sensitivity derived in 
Sect. 2.2, it is possible to set an exclusion plot. In the assumption of time 
modulation of Eq. (4), the experimental data allow to rule out (see Fig. 6) 
variations of the solar neutrino flux with low frequency and large amplitude. 
As expected, there is a rapid loss of sensitivity for u > 
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3 Maximum likelihood analysis 



An other option is to apply the maximum likelihood method to the time list 
of the ^^Ge candidate events in the solar runs, including (on the contrary of 
what happens in the standard Gallex/GNO analysis) a time-modulated term 
in the likelihood function. The general maximum likelihood method described 
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Fig. 6. 68% and 95% CL exclusion plots in the frequency/amplitude plane derived 
from the Lomb-Scargle analysis of GALLEX/GNO data 

in Ref. [17] can then be applied to the case of a production rate varying as 
in Eq. (4). In this condition the mean number of ^^Ge atoms produced in the 
target depends on the times of start and end of exposure (and not only from 
the total exposure time) and can be calculated as 



tsOE 



N, 



71 



n(t; no, A, i^, (^) e-^"^*'^"^-*) dt 



(5) 
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(see Eq. 12 of Ref. [17]), where A71 is the decay rate of ^^Ge. It is hence pos- 
sible to estimate from the solar run data the best-fit values for the four free 
parameters no. A, v and </) of Eq. (4). 

A known time modulation is actually expected in the solar neutrino data 
because of the seasonal variation of the Sun-Earth distance. Since the excen- 
tricity of the Earth's orbit is e=0.0167 and the perihelion is on January 4*^*, 
the expected parameters for the l/rf^ geometrical modulation in the solar neu- 
trino flux are: A ^ 2e = 0.033, v = \ y~^ and </?= 3.81 rad® . Although, at 
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least at the present stage, the ^^Ga data have not enough sensitivity to reveal 
such a low-amplitude modulation {A ■ hq ~ 2.3 SNU, to be compared with 
the overall Gallex/GNO error, ~ 5.5 SNU), this is anyway explored using the 
maximum likelihood method: the frequency and phase are fixed to the known 
values reported above while the amplitude is left as a free parameter. The 
best-fit coming from the combined analysis of the 123 Gallex/GNO solar runs 
is no=69.3±5.9 SNU and A=0.01±0.11. The amplitude is well consistent with 
what expected, but it is nevertheless consistent with zero, confirming that the 
presently reached sensitivity is not good enough to test this modulation^ . 
The data analysis was then repeated including two modulation terms in the 
likelihood: one, whose parameters are fixed, to account for the annual modula- 
tion and the other, whose parameters are left free, for other possible unknown 
periodicities. In this case the best-fit point for the 123 solar runs is (statisti- 
cal errors only) no=69.1±6.1 SNU, A=0.57±0.38, z/=13.9±0.4 y-^ (in good 
agreement with the frequency 13.59 y~^ reported in Ref. [5]) and ip= -0.3±1.0 
rad. The hypothesis of modulated solar neutrino fiux can be compared with 
the usual null hypothesis of stationary fiux (with the geometrical correction) 
using the likelihood ratio test [18]. If the constraint A = is applied (i.e. one 
makes the hypothesis of a constant neutrino fiux) the best-fit value is of course 
no = 69.3 ±4.1 SNU (statistical error only). If 

A = ^"^";^^ = °^ (6) 



is the likelihood ratio, the statistics —2 In A has the same distribution of x^(3) 
[18]; in the present case, —2 In A = 3.40 and the corresponding p- value is 33.4%. 
The experimental data are hence statistically consistent with the hypothesis 
of a time-constant neutrino interaction rate, in agreement with the results of 
the Lomb-Scargle analysis. 

In Fig. 7 it is shown the likelihood spectrum —2 In A obtained fixing the fre- 
quency and fitting all the other parameters from the data. It can be seen 
from the plot that the absolute maximum of the likelihood is not much larger 
than the other local maxima: therefore, no clear modulation frequencies can 
be derived from the data. Using the same likelihood ratio test, it is pos- 
sible to compare the hypothesis of stationary rate with the hypothesis of a 
time-modulated rate with frequency 13.59 y^^ [5]: in this case — 21nA=3.28 
(2 d.o.f.) and the corresponding p-value is 19.4%. 



^ The exposure needed to reduce the statistical error to A • no -^ 2.3 SNU would be 
of ~ 700 ton-y, to be compared with 270 ton-y of Gallex and GNO. Moreover the 
present systematic error, which is 2.5 SNU for GNO, should be further reduced in 
order to reach the required level of sensitivity. 



^ 9 






7 t 

6 

5 

4 F 

3 

2 

1 



_1 r I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
2.5 5 7.5 10 12,5 15 17,5 20 22.5 25 

Frequency (l /y) 

Fig. 7. Likelihood spectrum —2lnX{i/) from the Gahex/GNO dataset vs. frequency. 
All the other parameters are left free in the fit. The black dot marks the best-fit 
point. 

4 Conclusions 



The experimental data of the Gallex and GNO solar neutrino experiments 
(123 solar runs, for a total exposure time of 3307 days) have been searched 
for possible time modulations in the capture rate using the Lomb-Scargle pe- 
riodogram and the maximum likelihood methods. For the latter analysis, the 
time list of candidate '^^Ge events in the single runs has been taken into ac- 
count. In both cases no statistical evidence of time variations has been found: 
given the present sensitivity, the results are fully consistent with the expec- 
tation of a constant solar neutrino capture rate over the data taking period 
(> 10 y). Though this fact does not automatically exclude other hypotheses 
that predict a time-dependent rate, modulations with low frequency and large 
amplitude can be ruled out. 
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